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We show that in balhstic mesoscopic SNS junctions the period of critical current vs. magnetic 
flux dependence (magnetic interference pattern), /c(*&), changes continuously and non-monotonically 
from "I>o to 2$o as the length-to-width ratio of the junction grows, or temperature drops. In diffusive 
mesoscopic junctions the change is even more drastic, with the first zero of /c($) appearing at 3<l?o. 

The effect is a manifestation of nonlocal relation between the supercurrent density and superfiuid 
velocity in the normal part of the system, with the characteristic scale = hvp /i-KksT (ballistic 
limit) or = \/ hD /2-KkBT (diffusive limit), the normal metal coherence length, and arises due to 
restriction of the quasiparticle phase space near the lateral boundaries of the junction. It explains 
the 2$o-periodicity recently observed by Heida et al. (Phys. Rev. B 57, R5618 (1998)). We 
obtained explicit analytical expressions for the magnetic interference pattern for a junction with an 
arbitrary length-to-width ratio. Experiments are proposed to directly observe the $o 2$o- and 
$0 — > 3$o-transitions. 

It is well established that the electrodynamics of the superconductors is nonlocal on the scale of <^o , and so is 
the current-phase relation (since vector potential and superconducting phase both enter through one gauge-invariant 
combination, superfiuid velocity). In the normal layer of an SNS system the place of is taken by the normal metal 
coherence length, = hvF/2'KkBT (in the ballistic limit, L <C h, where li is the impurity scattering length), or 
= \/ D /2'n:kBT , (in the diffusive limit, U <C L) where D = vph/'i is the diffusion constant of electrons in the dirty 
metal . This length can obviously greatly exceed at low enough temperatures, and is limited only by the inelastic 
scattering length in the system, Z^, which also diverges as T ^ 0. (The latter gives the scale of classical nonlocality 
in mesoscopic systems, related to formation of essentially nonequilibrium states, when currents are determined by 
potential drops far from the observation point, like in classical point contacts or Landauer wires.) In wide SNS 
junction, the nonlocality is responsible for the formation of current-carrying Andreev levels, formed in the normal 
part of the system by Andreev reflections of quasiparticles from the off-diagonal scattering potential (order parameter) 
of the superconducting banks (^-|^, which provide the mechanism for the Josephson current flow. In mesoscopic SNS 
systems it is revealed in such effects as e.g. conductance oscillations in Andreev interferometers(Ref. |^ and references 
therein). It came nevertheless as a surprise, when Heida et al. ||^ reported that the critical current vs. magnetic 
flux dependence (magnetic interference pattern) in mesoscopic junctions, formed by Nb electrodes in contact with 2D 
electron gas of InAs heterostructure (S-2DEG-S junction), has a period 2<I>o — hc/e, instead of the standard $0 [!)• 
The authors of |0 posed a question whether this unexpected result is due to a nonlocal superconducting current-phase 
relation specific for their narrow (length-to-width ratio L/W ^ 1, instead of usual L/W <C 1), ballistic, mesoscopic 
junction. They also asked whether the diffusive character of scattering at NS boundaries (which was prominent in the 
experimental samples) is important for establishing such nonlocality. A heuristic formula proposed in |^ successfully 
reproduced the periodicity, but not the shape of the observed magnetic interference pattern. 

Since the nonlocality under question is already present in wide SNS junctions, where the magnetic interference 
pattern is nevertheless $o-periodic, the problem clearly requires a more thorough consideration. For example, the 
dephasing and inelastic scattering processes cannot be responsible for the "wide- narrow" difference, because the 
standard theory predicts <I>o-periodicity even when 1^,S,t = 00. 

In this paper we show that $0 2<I>o-transition is a size effect, resulting from the restrictions on allowed quasi- 
classical trajectories of quasiparticles in the normal layer near its lateral edges. It is present already in the absence of 
diffusive Andreev scattering (perfect NS boundaries). As the L/W ratio drops, relative contributions from the edges 
diminish, and in the limit L/W the standard ^o-periodicity is restored. The latter is actually a result of certain 
contributions to the current cancelling each other, not of their decay because of junction size exceeding the inelastic 
scattering length. We also predict even more drastic effect in the diffusive SNS junction, that can be qualitatively 
explained along the same lines. 

An important consequence of the "geometric" origin of the effect is the continuous transition of the periodicity 
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between $o and 2$o as the L/W ratio changes. Indeed, the restriction on contributing quasiparticle trajectories can 
be considered as having a smaller effective flux through the system. The effective flux, and therefore the observed 
periodicity, obviously will be a continuous function of L/W . This transition could be directly observed in a modifi- 
cation of the experimental setup of , where the width of the channel is determined by the voltage applied to gate 
electrodes (Fig.^). Such split-gate technique is now widely used in hybrid systems 2DEG-superconductor (see e.g. 
pot ) . Another experimental possibility, applicable to the junctions of fixed dimensions (like in ) , will be discussed 
later. 

First consider a ballistic S-2DEG-S junction (Figjl^) in the limit i assuming perfect Andreev refiections 

at NS interfaces, and totally absorbing lateral boundaries (at y — ±W/2). The latter condition, strictly speaking, 
corresponds to an "open" SNS junction (Fig.^), analogous to the structures investigated experimentally in In 
case of laterally restricted 2DEG one should use, e.g., condition of mirror reflection from the side walls. We will see 
that this does not qualitatively change the results. 

The magnetic field and the vector potential (in the Landau gauge) are 

H = i7i, A = Hyi. (1) 

Finally, we use the standard steplike approximation for the superconducting order parameter in the banks of the 
junction: 

A(r) = |Ao|e*^i0(-L/2 - z) + |Ao|e*'*^0(z - 1/2). (2) 

(In the gauge (^ the superconducting phase is not affected by the field.) 

The Fermi wavelength in the normal part of the system, Xp = h/pp (w 0.1 /im), is small compared to all other 
length scales in the problem. This allows us to apply the method of quasiclassical Green's functions integrated over 
energies p2| , 

5^(n,r) = -/ (^g^{^,n;r); /^(n,r)== / (^T^{^,n;r); 

/Ll(n,r) = /c^(-n,r)*; 
g^{n,rf + /^(n, r)/^(n, r)+ = 1. 

Here f/i^(^, n; r), J\j(^, n; r) are spatial Fourier transforms of normal and anomalous Gor'kov functions of arguments 
Xi_2 = r ± (5r. The momentum canonically conjugated to is 



p = \/2m*^n. 

The supercurrent density is then found as a sum over Matsubara frequencies, 

j(r) = -WAA(0)«fT^ (ng^(n,r)) = -4^zeAA(0)z;FT ^ (n(.g^(n, r) - .g^(-n, r)))„^>o ; (3) 

UJ>0 UJ>0 

the angular brackets denote angular averaging over the Fermi surface. 

Functions g^^, f^^ satisfy the set of Eilenberger equations, that in the ballistic limit and in the absence of the field is 
11 

(2t^ + UFn- V)^(n,r) = 2A(r).g„(n,r); (4) 
vpn ■ V5.(n, r) = A(r)*^(n, r) - A(r)/+(n, r). (5) 

Its characteristics are straight lines, interpreted as quasiclassical trajectories of quasiparticles It should be solved 
in three regions: z < —L/2, \z\ < L/2, z > L/2, using the limiting values in the bulk of left, right superconductor, 

Un,z^T^)= ^ ; (6) 

g^(n,z = too) = ^==^==, (7) 

and then matching solutions over the interfaces. Eventually, this will yield the standard sawtooth expression for the 
superconducting current density in an SNS junction which we will write in the form [ p^ for a point at the 

right NS boundary: 



2 



, ~ /" „efi? COS0 2 , ,1,,, L sin — 



Here ^^,(0) = Z!^^^ = cos 0. 

In our assumptions (totally absorbing walls) the integration over 9 must be limited to the directions within the 
angle shown in Fig.|l|, and the total Josephson current is written as 

j(, , N e^F [ [^'^ dyidy2 2^ L sinfc(0i-02) 

M^A,iil^/2[l + (^)2]3/^^^;" ^ M^..-.Jsinh,^^ 



(we take into account that tan By-^-y^ = {y2 — yi)/L). 

The only difference that nonzero magnetic field will make is that (0i — (^2) in (||) will be replaced by 



27r f° 

(p{n, r) ^ 4)2 ~ 4>i + — / A(r - vprn) ■ ndr, 
^0 Jt, 



(10) 



where the vector potential is integrated along the quasiclassical trajectory. (We neglect here the dynamical effects of 
the magnetic field, small by the parameter hu>c/ fJ--) The expression for the Josephson current becomes 

lU) = / 2^ L ^^^k{^iy,+y2) + <p) _ 

/c = max 1(6). (12) 

The current is given by a sum of the contributions from quasiclassical " Andreev tubes" of width ^ Xp each, following 
the quasiparticle trajectories. Each contribution depends on the length of the trajectory, and on the phase gained 
along it (both from Andreev reflections at the NS boundaries, and from the vector potential on the way through 
the normal part of the system). You will also notice that heuristic formula (7) of Ref. correctly captures the 
qualitative picture of the effect and follows from our expression in the limit of narrow junction at high temperature 
(L/W^,L/er^oo). 

At zero temperature, the results can be obtained explicitly in the limiting cases of wide (L/W 0) and narrow 
{L/W — > 00) junctions. Introducing 1/ = ^/^o, we obtain (see Figj^) 

«.,.?5:£Z£(i^MM, v,,,-^0; ,13, 

where the first formula is the standard result for a wide SNS junction; all harmonics contribute to the current, 
producing the distinctly non-Fraunhofer picture By {x} we denote the fractional part of a;. 

In general case, the dependence /c(i^) can be calculated numerically from (|l^ ) (see Fig.||). The period is indeed 
smoothly evolving from $0 to 2$o- Near the wide limit we see almost uniform growth of the period, in agreement 
with our qualitative considerations, but in the intermediate regime, L ^ W, the behaviour of /c(i^) becomes more 
complicated. It can be better understood in the high temperature limit, when only the term with fc = 1 in (|Tl| ) 
survives. Then the answer can be obtained explicitly for an arbitrary value of = L/W. Denoting by Dx the 
" effective aperture" , 

Dt = (15) 



we can write the expression for the critical current as 



^2 n2 „2 



fiu) = ly-^e^^^^^lm (Ei-i{{l/DT + iTTDTi^)/\/2) -Erf{iTrDTiy/V2)^ . (17) 
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The function /(i/) is proportional to sin(7ri/)/i/ at L/W = 0, and becomes strictly positive (proportional to 
sin^(7rj//2)/j/^) as L/W oo (Fig.^). Thus the behaviour of /c(i^) near its zeros shown in FigJ^ is due to merg- 
ing of real roots of /(i^). We see from (^,^^ that at finite temperature the "wide-narrow" transition is governed by 
mm{DT, L/W). Therefore another way of observation of the effect is by changing the temperature of the sample at 
fixed L, W. It has an advantage of being applicable to the mesoscopic SNS junctions of fixed dimensions. 

What difference would make a different boundary condition in the lateral direction? In case of mirror reflection from 
the side walls the current is given by a simple generalization of (jTl]), including the contributions from the reflected 
trajectories. In the absence of the magnetic field it is evident after we unfold the reflected trajectories by periodically 
extending the system in y-direction and taking the integral over yi from — oo to oo, that the Josephson current per 
width W of the contact will be the same as in an infinitely wide junction. But in the presence of the field the situation 
changes. Instead of linearly growing as Hy (Eq.(|^)), the effective vector potential in the extended system will be 
periodic in y ||T^, and the field contribution to the phase gain along the reflected trajectories will be systematically 
less than along their counterparts in an infinitely wide junction. Therefore the qualitative picture of the evolution of 
the oscillation period remains valid, as it is clearly seen in Fig.|[ 

Effect of strong elastic scattering in the normal layer is more drastic. Let us consider the limiting case of dirty 
SNS junction (FigJ^). Now the system is described by the Usadel equations for the quasiclassical Green's functions 
averaged over directions ||l7[| , F^^ and G^j- These equations and the expression for the supercurrent density read 
[Q(from now on we put h = c ~ 1): 

\u;\FUr) + ^[F^ir)V^GUr)-GUr){^ + 2ieA{r)fFUr)) = A*(r); (18) 

\F^\^ + Gl^l; (19) 
j(r) =-^eAA(0)OT^Im(^^*(r)(V + 2ieA(r))F^(r)). (20) 

We will limit our considerations to long enough junctions and high enough temperatures, L > ^t- Then the sum 
over Matsubara frequencies in (|o|) will be dominated by the terms with |a;| — nksT; the anomalous Usadel function 
will be small inside the normal layer. Therefore the first Usadel equation can be linearized to yield 

q^F^{r) - (V + 2ieA{r)f F„(r) = 0, (21) 

where = 2\uj\/D = 1/fy. The boundary conditions for F^{r) at z = ±L/2 and y = ±.W can be chosen as Q 

FUy,-L/2)^fo{cuy, FUy,L/2)^fo{cj)e'^^^-^'\ (22) 



where /o(w) — A±l/2/ {\uj\ + ■y/|wp + A±^^) . We are not interested in the actual value of the order parameter on the 
NS boundary (which due to proximity effect is lower than in the bulk) , as long as the anomalous Green's function stays 
small. We choose zero boundary conditions aX y = zLW/2, which would correspond to e.g. contact with clean normal 
conductor; the Josephson current is again carried only by the quasiparticles whose trajectories link the superconductors 
(Fig.^). After eliminating the vector potential from ( |20| , pl| ) by the transformation Fi^{r) — F^(r)e~^''^'^'", boundary 
conditions at zLL/2 are modified: 

FUy, -L/2) - h{Lo)e-'^-y; F^y, L/2) = /o(c.)e'(*^-^i)e*-''^. (23) 

The solution is readily expressed through the Green's function of Eq.(pl]), Q{y,z]rj,C,)^ in the rectangle with zero 
boundary conditions |1^ : 



27r .Lw/2 \ dz 



-e^(^--^-)e"^y ^g{y,z;fj,C) 



z=L/2 I 



The latter can be found by the method of mirror images: 

oo oo 



m— — oo n— — oo 



Here Q^{y,z;riX) — j^Ko{qT\/ {y — t\Y + ^ C)^) is the Green's function in the infinite plane, Kq is the modified 
Bessel function. 
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Summation over m can be limited to m = 0, ±1 because of our assumption L > ^t- The sum over n is taken using 
the Poisson summation formula. Substituting this in the expression for the current density at z = 0, we finally find 
for the Josephson current: 

I{(f>2 - 01, J^) = /c(j^) sin((/)2 - (pi), 

where 

= neAf{Q)DT\M^ , (25) 



Si{u) = VW^iqT + nHVw')'/'Ki/2{^u^iql + TT^iyW^)); SKu) = ^Si{u). (27) 

Behaviour of the function fdiff{i^) is shown in Fig.^. Due to exponential decay of Si{L/2) and its derivative, in the 
limit of narrow junction (W ^ ^t, L) only the terms with smallest / = ±1 survive This yields fdiff{t^)/fdiff{^) = 
cos^(7rz//2)/(l — v^Y . The first zero of this function is at $ = 3<I>o, tripling the first oscillation period. The following 
zeros are separated by 2$o-intervals (though the amplitude of higher-order peaks decreases so fast with that the 
possibility to observe this periodicity is rather problematic). In the opposite limit of wide junctions, W ':S> £,Tt Si's 
cease to depend on Z, and the shape of the interference pattern will be given by 



y / sanT{v + l)/2 _ ; sin7r(^ - /)/2 V 
^ ^ ' \ tt{v + 1)/2 ^ ' n{v^l)/2 ) 



1— — 00 



-4 , 

TTZ/ 



the Fraunhofer pattern expected in a wide dirty SNS junction. As in the ballistic case, continious change of the 
periodicity can be caused by lowering the temperature of the system. 

Note that in the range of parameters {W ^ L ^ ^t) when the first zero of the critical current appears at v ^ 2 the 
curve Ic{v) (oc \fdiff{i^)\) is distinctly sharper than in the ballistic case (see Figs.^^. Since the measurements of Ref. 
[0 (with diffusive scattering at NS boundaries) produced sharp cusps rather than smooth minima at $ ^ 2<I'o7 one 
can speculate that their system was effectively in diffusive rather than ballistic regime. 

Qualitatively the effect can also be understood in terms of loss of available phase space for current-carrying quasi- 
particles. Only the diffusive trajectories linking both superconductors contribute to the Josephson current (carried 
along such a trajectory by electrons and holes related through Andreev reflections at the end points). Therefore the 
trajectories originating too close to the lateral boundaries have more chances to end on the boundary and be lost, 
and the effective width of the junction, as well as the effective magnetic flux through it, are less than their actual 
values. Finite width of the junction also limits the "wandering" trajectories with large effective phase gain. Due to 
strong elastic scattering, the average z-component of the current carried through such a trajectory is the same as for 
a "shortcut" one, quite unlike the ballistic case where contributions from "grazing" trajectories are systematically 
less. This explains why the periodicity change is even more drastic in the diffusive case. 

In conclusion, we have demonstrated that the periodicity of the critical current dependence on the magnetic flux 
penetrating a mesoscopic SNS junction generally depends on the geometry of the system and is not given by either of 
fundamental quanta, $o = hc/2e or 2$o = hc/e, in a contrast to the case of tunneling Josephson junction. The effect 
is a manifestation of the nonlocal electrodynamics in hybrid normal-superconducting systems on scale Ct(Ct) S> Co- It 
can be understood in terms of quasiclassical Andreev levels ("Andreev tubes"), following the quasiparticle trajectories 
in the normal part of the system. Our theory provides an explanation for the recent experimental observation of 2<I>o- 
periodic magnetic interference pattern in S-2DEG-S junctions Q], and predicts a continuous transition from <I>o to 
2$o periodicity in ballistic case, and to 2$o-periodicity with the first iscillation period 3<I>o on the dirty limit. This 
transition can be observed either by changing the width of the 2DEG layer in the split-gate technique, or by changing 
the temperature of the fixed-size junction. The latter approach can be applied as well to mesoscopic SNS junctions 
with normal metal layer instead of 2DEG. 
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FIG. 1. (a) Josephson current in a ballistic S-2DEG-S junction, (b) "Open" SNS junction. Only the quasiparticle trajectories 
linking both superconductors contribute to the Josephson current, (c) Suggested experiment for direct observation of the 
$0 — > 24'0-transition. The width of the system is changed continuously by the gate voltage, Vg. 
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FIG. 3. Magnetic interference pattern and positions of zeros of the critical current as a function of length-to-width ratio (at 
L = 5^t)- Note the nonmonotonic behaviour of even zeros, and the fact that exact 3>o- or 2$o- periodicity takes place only in 
the limiting cases of infinitely wide (narrow) junction. 
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FIG. 6. Josephson current in an SNS junction in diffusive limit. Trajectory A does contribute to tlie Josephson current, 
trajectory B does not. 
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